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Abstract 

We  establish  a  new  group-contribution  model  based  on  the  Nernst-Einstein  equation  in  which  the  diffusion  coefficient  is  derived  from  the 
modified  double-lattice  (MDL)  model  and  the  Debye-Hiickel  (DH)  theory.  The  model  includes  the  combinatorial  energy  contribution  that  is 
responsible  for  the  revised  Flory-Huggins  entropy  of  mixing,  the  van  der  Waals  energy  contribution  from  dispersion,  and  the  polar  force  and 
the  specific  energy  contribution  from  hydrogen  bonding. 

The  Nernst-Einstein  equation  takes  into  account  the  mobility  of  the  salt  and  the  motion  of  the  polymer  host.  To  describe  the  segmental 
motion  of  the  polymer  chain,  which  is  the  well  known  conduction  mechanism  for  solid  polymer  electrolyte  (SPE)  systems,  the  effective 
co-ordinated  unit  parameter  is  introduced. 

Our  results  show  that  good  agreement  is  obtained  upon  comparison  with  experimental  data  of  various  PEO  and  salt  systems  in  the  interested 
ranges. 

©  2005  Elsevier  B  .  V.  All  rights  reserved. 
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1.  Introduction 

Solid  polymer  electrolytes  (SPEs),  formed  by  dissolving 
salts  in  a  polymer  matrix,  have  long  received  extensive  atten¬ 
tion  because  of  their  potential  for  achieving  high  ionic  con¬ 
ductivities.  Their  many  direct  applications  are  high-energy- 
density  batteries,  electrochromic  display  devices  (ECD),  and 
fuel  cells  [1-5]. 

Solid  polymer  electrolytes  have  improved  safety  for 
lithium  battery  compared  to  liquid  electrolytes,  but  they  are 
known  to  result  with  insufficient  performance,  especially  due 
to  low  ionic  conductivity.  Researches  on  the  electrochemical 
applications  of  SPEs  have,  therefore,  focused  on  the  ionic 
conductivity  for  each  complex. 
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To  increase  the  ionic  conductivities  of  SPE,  the  polymer 
should  have  both  low  Tg  and  low  crystallinity.  Reibel 
et  al.  [6]  described  ionic  conductivities  of  PEO/lithium 
bis(4-nitrophenylsulfonyl)imide  (LiNPSI),  and  PEO/lithium 
bis(trifluoromethanesulfonyl)imide  (LiTFSI)  systems  as 
polymer  electrolytes.  Sreekanth  et  al.  [7]  reported  an 
application  of  PEO  complexes  with  NaN03  salt  as  an 
electrochemical  cell.  Andreev  and  Bruce  [8]  investigated 
the  structure-conductivity  relation  of  PEO/LiAsF6  in  the 
analogous  phase.  Reddy  and  Chu  [9]  also  reported  the 
structure-conductivity  relation  of  PEO  with  potassium 
ionic  salt.  Historically,  for  the  quantitative  study  of  the 
ionic  conduction,  Gibbs  and  co-workers  [10,11]  developed 
the  configurational  entropy  model  for  polymer  properties, 
which  has  been  discussed  extensively  by  Goldstein  [12]. 
Configurational  entropy  (or  free  volume)  models  can  be 
expressed  in  the  form  of  a  Vogel-Tamman-Fulchur  (VTF) 
equation  [13-15],  which  is  generally  used  for  describing 
the  dependence  of  conductivity  on  temperature.  S0rensen 


J.H.  Joo,  Y.C.  Bae  /  Journal  of  Power  Sources  157  (2006)  448-456 


449 


and  Jacobsen  [16]  developed  a  simple  model  that  accounts 
quantitatively  for  the  concentration  dependence  of  the 
conductivity  for  the  low  purity  type  electrolyte.  MacFarlane 
et  al.  [17]  investigated  the  effect  of  plasticizer  on  the 
conductivity  by  using  the  Adam-Gibbs  model  combined 
with  Flory’s  configurational  entropy  at  the  fixed  temperature. 

Recently,  Kim  and  Bae  [18,19]  developed  configurational 
entropy  model  for  conductivities  of  SPEs  that  can  express 
the  composition  dependence  of  the  given  systems  based 
both  on  the  Adam  and  Gibbs  conductivity  model  and  on 
the  Flory’s  entropy  model.  To  take  into  account  the  pres¬ 
sure  effect  on  the  ionic  conductivities  of  the  compressed 
SPE  systems,  Ahn  and  Bae  [20,21]  and  Choi  and  Bae  [22] 
extended  Kim’s  configurational  entropy  model.  It  is  impor¬ 
tant  that  the  salt  effect  must  be  taken  into  account  in  the  ionic 
conductivity  model  because  the  ionic  conduction  is  worked 
by  the  mobility  of  salt  as  well  as  the  segmental  motion  of 
polymer. 

The  purpose  of  this  work  is  to  overcome  the  shortcomings 
of  the  previous  models  by  taking  into  account  the  salt  effect  in 
the  conductivity  model  using  a  group-contribution  method. 

Group-contribution  method  is  very  efficient  tool  to 
describe  thermodynamic  properties  of  polymer  solution 
because  it  utilizes  exiting  phase  equilibrium  data  when  pre¬ 
dicting  phase  behaviors  of  given  systems  of  which  data  are 
not  plentiful.  The  ultimate  goal  of  the  group-contribution 
method  lies  in  its  ability  to  predict  physical  properties  for 
systems  which  are  not  included  in  the  experimental  data,  that 
is  the  set  of  data  uses  to  determine  the  parameters.  The  basic 
idea  is  starting  from  that  whereas  the  chemical  compounds 
of  interest  in  chemical  technology  are  numerous,  the  number 
of  functional  groups  which  constitute  those  compounds  is, 
however,  much  smaller.  Thermodynamic  properties  of  a  fluid 
then  can  be  calculated  as  sum  of  contributions  made  by  the 
functional  groups.  However,  any  group-contribution  method 
is  necessarily  approximate  as  the  contribution  of  given  group 
in  one  molecule  is  not  necessarily  the  same  as  that  in  another 
molecule. 

To  describe  the  ionic  conductivity  of  the  given  systems,  we 
employ  the  diffusion  coefficient  for  which  the  driving  force 
is  based  on  a  gradient  of  chemical  potential. 

To  account  for  the  effect  of  both  salt  and  polymer  concur¬ 
rently,  the  sum  of  each  chemical  potential  is  differentiated 
with  concentration.  Since  the  conduction  mechanism  in  SPE 
is  correlated  with  segmental  motion  of  polymer  chain,  the 
effective  co-ordination  between  salt  and  polymer  units  is 
taken  into  account. 

Each  chemical  potential  is  calculated  from  Flory’s 
melting  point  depression  theory  [23]  combined  with 
the  modified  double-lattice  (MDL)  theory  and  the 
Debye-Hiickel  (DH)  theory  extended  by  Guggenheim.  Com¬ 
bination  of  the  derived  diffusion  coefficient  equation  with 
Nemst-Einstein  relationship  yields  the  final  conductivity 
equation. 

In  this  work,  we  develop  a  new  group-contribution  model 
taking  into  account  interactions  between  different  species 


(salt).  It  describes  ionic  conductivities  of  number  of  S PE/Li 
and  sodium  salt  systems  at  various  temperatures  and  compo¬ 
sitions. 


2.  Model  description 


Three  theoretical  aspects  are  taken  into  account:  the  lat¬ 
tice  notation  of  a  Debye-Hiickel  type  function  proposed  by 
Guggenheim  [25,26],  modified  double-lattice  model  [24]  and 
Flory’s  melting  point  depression  concept  [23]. 

In  this  study,  the  expression  for  the  Helmholtz  energy  of 
mixing  for  binary  polymer  solutions  is  defined  as  a  sum  of 
two  contributions. 

The  total  molar  Gibbs  energy  of  mixing  AG™  is 
assumed  to  consist  of  an  extended  Debye-Hiickel  theory  and 
modified  double-lattice  theory: 
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where  AG^L/RT  is  the  Helmholtz  energy  of  modified 
double-lattice  model;  A G^JRT,  the  Helmholtz  energy  of 
Debye-Hiickel  model;  R ,  the  gas  constant  and  T  is  the  abso¬ 
lute  temperature. 


2.7.  Modified  double-lattice  model 


2.1.1.  Primary  lattice 

Oh  and  Bae  [24]  proposed  a  new  Helmholtz  energy  of 
mixing  as  the  form  of  Flory-Huggins  theory.  The  expression 
is  given  by: 

Wr  =  (tt)  +  (|) ln*  +  (2) 


where  Nv  is  the  total  number  of  lattice  sites;  k,  the  Boltzmann’s 
constant  and  r*  is  the  number  of  segments  per  molecule  i.  xob 
is  a  new  interaction  parameter  and  functions  of  r,  and  e  and 
is  given  by: 


XOB  =  Cp 
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£  is  a  reduced  interaction  parameter  given  by: 

£  (£ii  +  £22  —  2£i2> 

£  =  —  =  -  (4) 

kT  kT 

where  £n,  £22  and  £12  are  for  the  corresponding  nearest- 
neighbor  segment-segment  interactions.  Parameters  Cp  and 
Cy  are  universal  constants.  These  constants  are  not  adjustable 
parameters  and  are  determined  by  comparing  with  Madden 
et  al.’s  Monte-Carlo  simulation  data  ( r\  =  1  and  r2  =  100). 
The  best  fitting  values  of  Cp  and  Cy  are  0.1415  and  1.7985, 
respectively  [24]. 
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2.1.2.  Secondary  lattice 

In  Freed’s  theory  [27,28],  the  solution  of  the  Helmholtz 
energy  of  mixing  for  the  Ising  model  is  given  by: 


AA 

NrkT 


=  x\  \nx\  +  X2  In X2  + 


ZSX\X2  Z£2x\x 2 


+ 


(5) 


where  z  is  the  coordination  number  and  Xi  is  the  mole  fraction 
of  the  component  i. 

To  obtain  an  analytical  expression  for  the  secondary  lat¬ 
tice,  we  defined  a  new  Helmholtz  energy  of  mixing  as  the 
fractional  form  to  improve  the  mathematical  approximation 
defect  by  revising  Eq.  (5).  This  secondary  lattice  is  introduced 
as  a  perturbation  to  account  for  the  oriented  interaction.  The 
expression  is  given  by: 
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where  A AseCjy  is  the  Helmholtz  energy  of  mixing  of  the  sec¬ 
ondary  lattice  for  i-j  segment-segment  pair;  Ny,  the  number 
of  i-j  pairs,  the  reduced  energy  parameter  contributed  by 
the  oriented  interactions  and  ij  is  the  surface  fraction  permit¬ 
ting  oriented  interactions.  For  simplicity,  ij  is  arbitrarily  set 
to  0.3  as  Hu  and  co-workers  [29,30]  suggested.  Ca  is  also 
not  an  adjustable  parameter  and  is  determined  by  comparing 
with  Panagiotopolous  et  al.’s  Gibbs-ensemble  Monte-Carlo 
simulation  data  of  Ising  lattice.  The  best  fitting  value  of  Ca 
is  0.4880  [24]. 

2.1.3.  Incorporation  of  secondary  lattice  into  primary 
lattice 

To  incorporate  a  secondary  lattice,  we  replace  Sy 
by  £y  —  AAseCty/Ny  in  Eq.  (4).  If  oriented  interaction 
occurs  in  the  i-j  segment-segment  pairs,  we  replace  e  by 
£/kT+2(AASQCjj/NijkT)  in  Eq.  (5).  If  oriented  interaction 
occurs  in  the  i-i  segment-segment  pairs,  we  replace  e  by 
s/kT  —  A Asecji/NakT.  In  this  study,  we  assume  the  oriented 
interaction  occurs  in  the  i-i ,  j-j  and  i-j  segment-segment 
pairs.  We  replace  e  by: 
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where  s\v  s\2  and  s\2  are  van  der  Waals  energy  interaction 
parameters.  A AseCjn,  AAsec?22  and  AAsec  j2  are  the  additional 
Helmholtz  functions  for  the  corresponding  secondary  lattice. 
Eq.  (7)  then  becomes: 
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To  correlate  MDL  model  to  melting  point  depression  the¬ 
ory,  we  require  chemical  potentials  of  components  1  and  2. 
The  definition  of  chemical  potential  is: 
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The  final  expression  for  the  chemical  potential  can  be  writ¬ 
ten  as: 
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where  0/  is  the  segment  fraction  of  component  /,  0/  =  NirfNr 
and  Nr  =  YU? Niri  is  the  total  number  of  segments  in  the 
system  and  r;  is  the  segment  number  of  components  1  (salt) 
and  2  (polymer). 


+ 


$£22/ kT 
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2.2.  A  lattice  notation  of  an  extended  Debye-Huckel 
theory 


For  a  binary  polymer/salt  system  at  solute  molality  m 
(mol  kg-1  polymer),  Guggenheim’s  expression  for  the  molar 
Gibbs  energy  of  mixing  AG™ could  be  rewritten  in  the 
framework  of  lattice  theory  as  follows  [25,26]: 


AG™  01 


RT 


r\vm 


ln(l  v)  —  x  H- 


(12) 

(13) 


where  <fi\  is  the  segment  fraction  of  the  salt  ion,  r\  (=  1),  the 
number  of  segments  per  salt  ion,  v(=  vm  +  vx\  where  r>M 
and  vx  are  the  number  of  M  and  X  ions  per  salt,  respec¬ 
tively),  the  number  of  ions  per  salt  and  /  is  the  ionic  strength. 
A  is  the  usual  Debye-Huckel  coefficient.  In  this  study  we 
fix  A  =  0.068  (assumed  to  be  independent  of  temperature)  as 
an  optimization  factor  for  polymer/salt  systems.  This  small 
value  implies  that  ion-ion  interactions  in  a  polymer/salt  sys¬ 
tem  are  relatively  small  (for  example,  calculated  percentage 
at  molality  of  salt  ^0.1  mol  kg-1  in  PE0/LiCF3S03  system: 
ions  =  2%;  pairs  =  71%;  and  triples  =  27%)  [1].  For  a  binary 
polymer/salt  system  containing  1  kg  of  polymer  and  vm  moles 
of  salt  ion,  0 1,  02  and  I  are  defined  by: 
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where  M  is  the  molecular  weight  of  polymer  in  g  mol-1  (i.e. 
M-  900,000  g  mol-1).  Afi,  Afo,  zm,  zx  and  r2  are  the  number 
of  moles  of  salt  ion  and  polymer,  the  valences  of  M  and  X 
ions  and  the  number  of  segments  per  polymer,  respectively. 
The  chemical  potentials  are  given  by: 
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2.3.  van  der  Waals  energy  contribution 


The  energy  parameter  £*•  in  Eq.  ( 1 9)  is  due  to  van  der  Waals 
forces  (dispersion  and  polar  forces).  For  a  pure-component  i, 
£*•  can  be  estimated  using  the  square  of  the  pure-component 
van  der  Waals  solubility  parameter  of  Hansen  (Barton)  [31], 
which  is  the  sum  of  a  dispersion  contribution  and  a  polar 
contribution:  82dw  =  8d  +  82. 


3NAs*iri 


(19) 


where  NA  is  the  Avogadro  number  and  where  <52dw  and  Vml 
are  at  25  °C.  For  a  pure-component,  the  effect  of  temperature 
on  £*;  is  given  by: 


£ii  —  0si£l  +  0  s2£2 


(20) 
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where  Vmi  depends  on  temperature.  The  temperature- 
independent  parameter  s ^  can  be  estimated  by: 
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15.17  x  10  6  m3  mol  1 


(22) 
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The  constant  15.17  x  10~6m3mol-1  is  the  molar  hard¬ 
core  volumes  of  a  CH2  group.  In  our  model,  the  group- 
contribution  concept  is  considered  to  calculate  the  chain 
length  contrary  to  that  of  the  existing  modified  double-lattice 
model. 

The  cross  interaction  van  der  Waals  energy  parameter  £*■  is 
estimated  by  the  geometric  mean  of  the  corresponding  pure- 
component  parameters: 


(24) 


Cross  specific  energy  parameter  bsij  is  calculated  from 
pair-interaction  group  parameters: 

fmfngmn  (25) 

where  Ns  and  Afp  are  number  of  groups  in  solvents  and  poly¬ 
mers,  respectively.  0m  and  fn  are  volume  fractions  of  group 
m  in  a  solvent  and  that  of  group  n  in  a  polymer,  respectively; 
and  gmn  are  pair  interaction  parameters  between  group  m  in 
a  solvent  and  group  n  in  a  polymer.  To  improve  the  accuracy 
of  prediction,  we  assume  that  a  functional  group  in  a  poly¬ 
mer  is  different  from  that  in  a  solvent.  In  this  study,  they  are 
estimated  by  fitting  experimental  solid-liquid  equilibria  data 
of  polymer  solutions. 


m  —  1  n = 1 
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2.4.  Ionic  conductivity 


ordinating  polymer  units,  which  is  given  by: 


For  binary  diffusion  in  gases  or  liquids,  the  generalized 
Fick’s  equation  for  heat  and  mass  is  as  follows  [32]: 


4  =  ~cDAB 
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xAV  In  a  a  +  — ~[(4>a  ~  oja  Wp 

cRT 


-  pdJAmigk  ~  gB)]  +  &tV  In  T 


(26) 


This  equation  represents  that  the  thermodynamics  of  irre¬ 
versible  processes  dictates  using  the  activity  gradient  as  the 
driving  force  for  concentration  diffusion.  This  requires  a  dif¬ 
fusion  coefficient  different  from  Fick’s  first  law.  When  the 
pressure-,  thermal-,  and  forced-diffusion  terms  are  dropped, 
Eq.  (26)  for  binary  electrolyte  is  simplified  by: 


7*  =  -D*CSV  In  as  (27) 

where  D  \  Cs  and  as  are  self-diffusion  coefficient,  concen¬ 
tration  and  activity  of  salt,  respectively.  This  equation  may 
be  rewritten  by  making  use  of  the  fact  that  the  activity  is  a 
function  of  concentration  to  obtain: 


r.  _  d(A Heff/RT) 

D  s  —  Lf  L,  s  , 

dCs 

A^eff  Apts  ,  .  Aptu 

-  —  -  i  Agff -  (51) 

RT  RT  RT 

where  Aeff  are  the  effective  co-ordinated  units  of  polymer.  A 
mathematical  form  of  AoewCs  is  adopted  for  Aeff  based  on  the 
exponentially  lowering  coordinating  units,  where  Ao  and  co 
are  adjustable  model  parameters,  respectively. 

This  in  turn  yields,  using  the  Nernst-Einstein  relationship 
for  multi-component  system,  the  conductivity  equation  for 
SPE  having  the  form  as: 

Y 2 

*  =  (32) 

i 

where  F  is  a  Faraday  constant.  Since  we  assume  that  the 
phase  at  the  given  condition  is  binary  system  of  polymer  and 
salt,  the  moving  object  is  salt  itself  instead  of  cation.  If  the 
charge  effect  of  each  ion  remains  for  the  ion  interactions,  this 
assumption  brings  Eq.  (32)  in  a  simple  form: 


*  *  /  d  In  as  \ 

j;  =  -D*Cs\—^)  VCs  (28) 

when  comparing  Eq.  (28)  with  the  original  Fick’s 
equation,/*  =  —  DSVCS,  this  is  related  to  the  measured  dif¬ 
fusion  coefficient  Ds  (based  on  a  concentration  driving  force) 
by  [32]: 


d  In  a, 
d  Inc. 


(29) 


where  D  characterizes  the  component  mobility  in  the 
absence  of  any  interactions  in  the  given  system  [33].  This 
may  be  rewritten  by  the  fact  that  the  activity  is  related  to  the 
chemical  potential  by  In  a  =  A /jl/RT 


CKAPW^T) 

dCs 


(30) 


The  transport  of  cations  in  solvent-free  polymer  elec¬ 
trolytes  differs  from  that  of  systems  based  on  molecular 
liquids  or  low  molar  mass  polymers.  In  the  latter  systems,  ions 
can  move  together  with  their  co-ordinated  solvent,  but  in  the 
case  of  high  molecular  weight  polymers  the  centre  of  grav¬ 
ity  of  the  chain  cannot  be  moved  significant  distances.  For 
electrolytes  using  high  molecular  weight  polymers,  cation 
transport  does  not  only  occur  in  conjunction  with  polymer 
diffusion  but  lithium  ion  transport  also  occurs  as  a  result 
of  segmental  motion  [34].  Since  cations  move  across  co¬ 
ordinating  sites  which  are  made  up  of  the  acid-base  inter¬ 
actions  between  solvent  and  solute  molecules,  diffusion  of 
cation  must  be  taken  into  account  with  effective  co-ordinated 
polymer  units.  To  express  this  co-ordinating  effect  in  the  con¬ 
ductivity  model,  the  chemical  potential  in  Eq.  (30)  is  replaced 
by  the  sum  of  chemical  potentials  of  salt  and  effective  co- 


a  — 


F2C, 

RT 


Ds 


(33) 


Substituting  Eq.  (30)  into  Eq.  (33)  gives  the  final  ionic 
conductivity  equation  for  SPE  systems: 


a  — 


F2Cs  d(Ans/ RT  +  Xoe^i AhJRT))^  2 

D  C, - - - 2^' 


RT 


d  G 


(34) 


where  the  chemical  potentials  are  given  by: 


A/Xs 

RT 


v  ( A\zmzx\  Il/2\ 

1  + 11'1  J 


1000 


+  ln(l  -  02)  -  n  (  —  -  —  )  02 

n 


+  r\ 


C, 


1 

r2 


1 


+  I  2  — | - l  £ 

r2 


1 

n 


02 


1 

r2 

2 


+ 


1 

r2 


1 

n 


+  Cy£  )  £ 


—  2r\ 


1  1 

r2  r\ 


- -  )  +  Cy£  )  £  +  Cy£ 


y 


02 


+  3nCys24) 


(35) 


and 
A/xu 

~RT  '  V\  r2  V  1000 


vun  (  m 


^A/3/2a(/1/2)j  +  In  02 


+  r2 


---)  +  (2+-  I  £ 

J2  r\J  \r2  r\J  V  r2 
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+  2(2+  —  I  s 
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02 


+  T2 


+  c. 
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02 
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(36) 


3.  Results  and  discussion 

We  establish  a  new  group-contribution  model  based  on 
the  Nernst-Einstein  equation  in  which  the  diffusion  coeffi¬ 
cient  is  derived  from  the  modified  double-lattice  model  and 
the  Debye-Hiickel  theory  with  A  =  0.068.  The  new  model 
employs  the  secondary  lattice  concept  to  take  into  account 
an  oriented  interaction.  The  advantage  of  this  model  follows 
from  its  simplicity. 

In  this  work,  most  of  parameters  are  calculated  from 
pure-component  properties,  either  from  experimental  data  or 
from  published  estimation  methods.  To  establish  the  group- 
contribution  method,  the  most  significant  role  is  to  determine 
the  cross-pair  interaction  between  polymer  and  salt  segments. 

In  energy  parameters  of  polymer/salt  systems,  ions  can 
move  together  with  their  co-ordinated  solvent,  but  in  the  case 
of  high  molecular  weight  polymers  the  centre  of  gravity  of  the 
chain  cannot  be  moved  significant  distances.  For  electrolytes 
using  high  molecular  weight  polymers,  cation  transport  does 
not  only  occur  in  conjunction  with  polymer  diffusion  but 
lithium  ion  transport  also  occurs  as  a  result  of  segmental 
motion  [34].  Since  cations  move  across  co-ordinating  sites 
that  are  made  up  of  the  acid-base  interactions  between  sol¬ 
vent  and  solute  molecules,  diffusion  of  cation  must  be  taken 
into  account  with  effective  co-ordinated  polymer  units. 


Table  2 


Diffusion  and  co-ordinated  unit  parameters  for  PEO/salt  systems 


T  (K) 

D  ( cm2  s  1 ) 

A.0 

co  (cm3  mol  1 ) 

338 

2.85  x  10“7 

15.713 

-275.247 

PEO/LiAsF6 

358 

6.23  x  10“7 

18.527 

-212.359 

378 

1.24  x  10“6 

17.677 

-192.672 

398 

2.21  x  10“6 

13.163 

-157.232 

338 

5.3  x  lO"8 

12.612 

-387.425 

PE0/LiC104 

358 

2.5  x  lO"7 

9.252 

-326.247 

378 

4.0  x  10“7 

8.347 

-276.437 

398 

8.0  x  10“7 

6.403 

-245.623 

PEO/NACF3SO3 

358 

8.2  x  10"8 

18.527 

-175.241 

Fig.  1.  The  calculated  value  of  diffusion  coefficient  (Z+)  plotted  against  T. 
D*  =  —9.3126  x  10“6  + 2.7456  x  10~ST. 


Table  1  gives  physical  properties  of  each  component  such 
as  melting  temperature,  heat  of  fusion,  molecular  weight, 
density,  and  molar  volume  [35]. 

We  assume  that  Ao  and  co  are  co-ordinating  unit  parame¬ 
ters.  D  is  a  diffusion  coefficient  dependent  on  temperature 
and  is  obtained  from  elsewhere  [37].  These  parameters  (Ao, 
co  and  D  )  are  reported  for  the  co-ordinating  effect  at  vari¬ 
ous  temperatures  (Table  2).  Pai  and  Bae  have  [37]  set  these 
parameters  constant. 

We  employ  the  same  assumption  as  Pai  and  Bae  [37]. 
These  parameters  (Ao,  co  and  D  ),  however,  are  not  adjustable 
parameters  but  universal  parameters  dependant  on  tempera- 
ture.  In  Figs.  1-3,  calculated  values  of  Ao,  co  and  D  are 


Table  1 


Physical  properties  of  PEO  and  salts 


On  (K) 

AH  (J  mol  1 ) 

M.W.  (g  mol-1) 

Density  (gem  3) 

Vu  (cm3  mol  1 ) 

PEO 

338.15 

6798. 00a 

5000000 

1.21 

36.6 

338.15 

8284. 32a 

900000 

1.21 

36.6 

NaCF3S03 

527.15 

10433.72 

172.06 

1.13 

36.0 

LiC104 

509.15 

14600.00 

106.39 

2.43 

43.8 

LiAsF6 

525.63 

57188.20 

195.85 

2.65 

73.9 

a  Junir1. 
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Fig.  2.  The  calculated  value  of  co-ordinating  units  parameters  (Ao)  plotted 
against  T.  7q  =  —586.04  +  3.328x  —  0.00458x2. 

compared  with  experimentally  determined  values  calculated 
by  Pai  and  Bae.  To  simply  correlate  with  simulation  data,  Ao, 
co  and  D  are  suggested  as  follow: 


Fig.  4.  Ionic  conductivity  as  a  function  of  weight%  of  NaCF3S03  at 
358.15  K.  (A)  The  solid  line  is  calculated  from  D* ,  an  adjustable  model 
parameter  and  Ao  and  co  are  calculated  by  Eqs.  (38)  and  (39);  (B)  the  dashed 
line  is  calculated  with  fixed  D  at  3.5  x  10-7  and  Ao  and  co  are  calculated  by 
using  Eqs.  (38)  and  (39);  and  (C)  the  dotted  line  is  calculated  from  the  final 
ionic  conductivity  equation  (Eq.  (34))  using  Eqs.  (37)-(39)  for  the  given 
systems.  The  dark  circles  are  experimental  data  reported  by  Hu  et  al.  [30], 
and  the  lines  are  calculated  values  using  the  proposed  model. 


D*  =  -9.31  X  10"6  +  2.75  x  1(T87 
co  =  17988.91  -  98.877  +  0.1372 
A0  =  -586.04  +  3.337  -  0.0045872 


(37) 

(38) 

(39) 


salt  concentration  is  calculated  from  the  salt  mass  fraction  as 
follows: 


Pirn 

Mi 


(40) 


These  Eqs.  (37)-(39)  provide  a  simple  and  accurate  the 
expression  for  predicting  ionic  conductivity. 

Fig.  4  shows  ionic  conductivities  of  PE0/NaCF3S03  sys¬ 
tem.  The  dark  circles  are  experimental  data  reported  by  Ma 
et  al.  [36]  and  the  lines  are  calculated  by  the  proposed  model. 
To  differentiate  chemical  potential  with  concentration,  the 


where  w\  is  the  weight  fraction  of  salt.  We  set  r\-  0.3  and 
z  =  6  as  suggested  by  Hu  et  al.  [29].  The  density  of  salt,  pi, 
then  is  given  by: 


1 

Pi 


1 

— n  +  w\ 

Pi 


(41) 


Fig.  3.  The  calculated  value  of  co-ordinating  units  parameters  (co)  plotted 
against  T.  co=  17988.91  -  98.87T+0.133T2. 


where  p^  is  the  density  of  PEO. 

In  Fig.  4,  each  line  is  calculated  based  on  the  different 
cases  of  Ao,  co  and  D  parameters.  The  dotted  line  is  calcu¬ 
lated  from  the  final  ionic  conductivity  equation  (Eq.  (34)) 
using  Eqs.  (37)-(39)  for  the  given  systems.  Dashed  line  is 
calculated  with  fixed  D  at  3.5  x  10  .  Ao  and  co  are  cal¬ 

culated  by  using  Eqs.  (38)  and  (39),  respectively.  Results 
show  that  these  two  cases  slightly  overestimate  ionic  con¬ 
ductivities  of  the  given  system.  The  significance  of  such 
disagreement  is  difficult  to  assess  because  D *  especially  is 
very  sensitive  to  the  given  system  and  also  is  very  depen- 
dant  on  temperature.  When  we  set  D  is  an  adjustable  model 
parameter  and  Ao  and  co  are  calculated  by  Eqs.  (38)  and 
(39),  respectively,  we  obtain  the  best  fitting  line  (solid  line  in 
Fig.  4). 

Figs.  5  and  6  show  ionic  conductivities  of  PEO/LiClOq 
and  PEO/LiAsF6  systems  in  which  the  experimental  data  are 
reported  by  Robitallie  and  Fauteux  [38].  These  figures  are 
also  calculated  by  the  ionic  conductivity  equation  (Eq.  (34)) 
with  Eqs.  (37)-(39).  Dotted,  dashed  and  solid  lines  are  also 
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co 
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0.1  0.2 
Wt.  %  of  LiCIO, 


0.1 


0.01 


IE-4 


IE-5 


IE-6 


0.00  0.05  0.10 


0.15  0.20  0.25 

Wt.  %  Of  LiCI04 


0.30  0.35 


Fig.  5.  (a)  Ionic  conductivity  as  a  function  of  weight%  of  LiC104  at 
378.15  K.  (A)  The  solid  line  is  calculated  from  D* ,  an  adjustable  model 
parameter  and  ko  and  co  are  calculated  by  Eqs.  (38)  and  (39);  (B)  the  dashed 
line  is  calculated  with  fixed  D  at  3.5  x  10~7  and  Ao  and  co  are  calculated  by 
using  Eqs.  (38)  and  (39);  and  (C)  the  dotted  line  is  calculated  from  the  final 
ionic  conductivity  equation  (Eq.  (34))  using  Eqs.  (37)-(39)  for  the  given 
systems,  (b)  Ionic  conductivity  of  PEO/LiCICU  system.  Experimental  data 
are  reported  by  Barton  [31]  and  the  lines  are  calculated  by  the  proposed 
model. 


Wt.  %  of  LiAsF6 


0.00  0.05  0.10  0.15  0.20  0.25  0.30  0.35  0.40 

Wt.  %  of  LiAsF6 

Fig.  6.  (a)  Ionic  conductivity  as  a  function  of  weight%  of  LiAsF6  at 
378.15  K.  (A)  The  solid  line  is  calculated  from  D* ,  an  adjustable  model 
parameter  and  ko  and  co  are  calculated  by  Eqs.  (38)  and  (39);  (B)  the  dashed 
line  is  calculated  with  fixed  D  at  3.5  x  10-7  and  Ao  and  co  are  calculated  by 
using  Eqs.  (38)  and  (39);  and  (C)  the  dotted  line  is  calculated  from  the  final 
ionic  conductivity  equation  (Eq.  (34))  using  Eqs.  (37)-(39)  for  the  given 
systems,  (b)  Ionic  conductivity  of  PEO/LiAsF6  system.  Experimental  data 
are  reported  by  Barton  [31]  and  the  lines  are  calculated  by  the  proposed 
model. 


4.  Conclusion 


evaluated  by  the  same  procedure  as  in  Fig.  4  with  previously 
obtained  group  interaction  energy  parameters  (see  Table  3.) 
with  no  additional  group  interaction  parameters.  As  shown 
in  Figs.  5(b)  and  6(b),  the  fitted  values  agree  fairly  well  with 
the  experimental  data. 

Table  3 


Group  interaction  parameter  gmn  (k) 


Salts 

Polymer  (PEO) 

ch2 

O 

Li 

-386.47 

-762.45 

CF3SO3 

869.91 

240.719 

Na 

-3518.75 

-3483.71 

AsF6 

235.05 

-495.72 

CIO4 

1243.33 

-257.99 

We  establish  a  new  group-contribution  model  based  on 
the  Nemst-Einstein  equation  in  which  the  diffusion  coeffi¬ 
cient  is  derived  from  modified  double-lattice  model  and  the 
Debye-Huckel  theory  with  A  =  0.068.  The  proposed  model 
gives  a  simplified  and  improved  expression  for  the  ionic  con¬ 
ductivity  for  polymer/salt  systems  that  includes  the  van  der 
Waals  energy  contribution.  The  ultimate  goal  of  this  work  is 
to  describe  ionic  conductivity  of  SPE/salts  system  using  the 
group-contribution  method. 

When  D  is  determined  by  Eq.  (37),  calculated  values  are 
slightly  overestimated  experimentally  observed  ionic  con¬ 
ductivities  of  the  given  system.  The  significance  of  such 
disagreement  is  difficult  to  assess  because  D *  especially  is 
very  sensitive  to  given  system  and  also  is  very  dependant 
on  temperature.  When  we  set  D  be  an  adjustable  model 
parameter  and  Ao  and  co  are  calculated  by  Eqs.  (38)  and  (39), 
respectively,  the  best  fitting  results  are  obtained. 
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To  characterize  most  common  SPE/salt  system,  more 
experimental  data  are  required  to  obtain  numerous  group 
parameters  for  other  electrolyte  systems,  and  then  to  extend 
the  group-contribution  method  to  a  larger  variety  of  systems. 
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